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CDF recently reported an excess of events in the /i7 missing Et {^t) channel that disagrees with 
the Standard Model prediction. No such excess was observed in the e^flj, channel. We explain 
the excess via resonant smuon production with a single dominant R-parity violating coupling A211, 
in the context of models where the gravitino is the lightest supersymmetric particle. The slepton 
decays to the lightest neutralino and a muon followed by neutralino decaying to a gravitino and 
photon. We determine a viable region of parameter space that fits the kinematical distributions 
of the Run I excess and illustrate the effect by examining the best fit point in detail. We provide 
predictions for an excess in the and photon channel at Run I and Run If. Run II will decisively 
rule out or confirm our scenario. 



I. INTRODUCTION 



CDF has recently presented results on the production of combinations involving at least one photon and one lepton 
(e or /i) in pp collisions at ^/s — 1.8 TeV, using 86 pb~^ of Tevatron 1994-95 data ||l|. In general the results were 
consistent with the Standard Model (SM), however 16 photon-lcpton events with large were observed, with 7.6±0.7 
expected. Moreover, 11 of these events involved muons (with 4.2 ± 0.5 expected) and only 5 electrons (with 3.4 ± 
0.3 expected), suggestive of a lepton flavour violating asymmetry involving muons. 

What can such a process be? A natural framework with explicit flavour violating couplings is provided by R- violating 
supersymmetry Q , which contains operators with a complicated flavour structure in the superpotential 

Wrpv = ^XijkL^LjEk + XijkLtQjDk + ^X'^jkUiD^Dk + HiUH2 (1) 

where L [Q) are the left-handed lepton (quark) superfields while E, D, and U contain the corresponding right-handed 
fields, and i,j,k generation indices. The second of the above terms is of particular interest, since it can lead to 
resonant slepton production in hadron-hadron collisions ^ , via the diagram that appears below. 

Such a resonance would lead to enhanced cross sections with a rich final state topology, which, as we are going 
to show, can explain the CDF anomaly. What would then be the structure of the associated operator? R-violating 
couplings have upper bounds coming from various flavour-violating processes Q]. Therefore, to get the requisite 
number of events to explain the observed anomaly, a sizable cross section is required which would then imply a 
process with valence quarks in the initial state. Since the events are seen in the muon channel, the operator can 
be specified to be L2Q1D1, which generates the couplings flud and Vfj^dd (and charge conjugates), along with other 
supersymmetrised copies involving squarks. This coupling, A211, is constrained from R^^ = T{tt — > ei^)/(7r — * /ij/) Q| 
to be < 0.059 X |. 




FIG. 1: Resonant smuon production and subsequent decay 
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Upon production, sleptons (in our case, smuons or sneutrinos) can in general decay via a large variety of channels 
Q if they are kinematically accessible. However, the crucial observation is that R-violating supersymmetry by itself 
may not account for the observed anomaly, because of the fact that the anomaly is observed in a channel where a 
photon is produced. However, if the gravitino (present in all models where supersymmetry is gauged) is the lightest 
supersymmetric particle (LSP) it is too long-lived to decay within the detector Thus, the gravitino, G, provides 
the missing energy signature since it is electrically neutral and interacts rather weakly with matter. If the neutralino, 
as is often the case, is dominantly photino, then the decay Xi ~^ can dominate M. It is interesting to note that 
the 6677^3^ event recorded by CDF Q can be explained by such a decay j^. 

Since at the moment there is neither enhancement in the two-fermion final state, nor observation of chains of cascade 
decays, the most natural explanation is that the R-conserving decay mode of the smuon which produces the lightest 
neutralino dominates over the rest while subsequently Xi ~^ ^7. The competing R-parity violating decay modes of 
Xi ~^ ^ii ^-nd Xi ~^ leading to njj^x or /z/xjj final states become negligible (as is the case here) when X'211 and 
are both small enough. Smuon decay into two jets via the R-parity violating mode is essentially unobservable 
because of the huge 2 jet background. For example, for a resonance mass of 200 GeV, only a a.B > 1.3 x 10^ pb is 
excluded at 95% C.L. |jl^. This will not provide a restrictive bound upon our scenario. 

It is worth stressing the clarity of the signatures, but also of future predictions in the case of a resonant process. 
Moreover, the presence of both slepton and sneutrino resonances are in principle to be expected, and we provide a 
prediction for "f^x events. The higher statistics in Run II of the Tevatron should allow verification our model. 

The new aspect of the model we present here compared to previous studies of resonant slepton production at hadron 
colliders , is to marry the gravitino LSP scenario with R-parity violating supersymmetry. This marriage has been 
considered before in the context of dark matter 



II. MODEL AND RESULTS 

We use the ISASUSY part of the ISAJET7.58 package [|l2| to generate the spectrum, branching ratios and decays of 
the sparticles. For an example of parameters, we choose (in the notation used by ref. |l^) A211 = 0.01, = 10~^ 
eV, tan/3 = 10, At^r.b ~ 0, and scan over the bino mass Mi and the slepton mass mj = ^ = mg^ ^ GeV. The 
values of A211 and 7713/2 are dictated by the need to have the decays shown in Fig. |l| being dominant. However, there 
are ranges of values in the i?- violating coupling and the gravitino mass where this decay chain is obtained. In fact, 
the acceptable ranges are an order of magnitude in X'211 ^^'i orders of magnitude in m^/2- fJ- together with other 
flavour diagonal soft supersymmetry breaking parameters are set to be so heavy that any superparticles except the 
first two generation sleptons, the lightest neutralino and the gravitino are too heavy to be produced or to contribute 
to cascade decays in Tevatron data. They therefore do not appear in this analysis. We have checked that this is true 
over a large volume of parameter space. We emphasise that this is a representative hyperplane in the supersymmetric 
parameter space and not a special choice. 

We use HERWIG6 . 3 including parton showering (but not including isolation cuts) to calculate cross-sections for 
single slepton production. A 7-in-active-region cut requires that the photon not have rapidity \ri\ > 1 or \ri\ < 0.05. 
The region 0.77 < 77 < 1.0, 75° < (j) < 90° is also excluded because it is not instrumented. Fiducial photon detection 
efficiency was set to be 81%, whereas for the muons it is 66% for 1.0 > |?7^| > 0.6 and 45% for \r]\ < 0.6. and the 
Et of both the muon and photon were required to be greater than 25 GeV. 

We calculate the difference in log likelihood between our model and the SM given by each kincmatical variable 
that was presented in ref. |jl|. This provides 95% C.L. limits upon Mi and Am = — Mi. We show the viable 
regions for the energy distributions: Et{j, yu), -^t, the mass m and transverse mass Mt distributions m^^, Mt{h]^t)j 
Mt{^t1): Mxi'jfJ-^T) various transverse angular separations A(j)ij, where i,j = fi,j,^rp. AR, defined as the 
distance in 77 — (/) space between the muon and the photon, is also used. It is not possible to take correlations between 
these different kinematical variables into account because we do not possess the multi-dimensional data. Therefore 
we resort to examining each one in turn and see to what extent each region overlaps. Fig. || shows that all of the 95% 
confidence level regions overlap at Mi « 90 GeV, Am — 25 — 40 GeV, indicating that our model is in good agreement 
with all of the observed kinematical properties of the events. The region at the bottom the plots is ruled out by LEP2 
from neutralino pair production Q. 

The most discriminating kinematical variable is Et{ij): which favors our model over the SM at the 3.3(t level at 
the best fit point Mi = 87 GeV and Am = 35 GeV. We refer to this point as "the best fit point" from now on, 
and examine its properties more closely. We show the predicted distribution of lepton Et in Fig. ^ and compare it 
with the excess of the data over the SM background. Important features of the sparticle spectrum are displayed in 
Table |. We also show the range of sparticle masses corresponding to the acceptable fit range of parameter space. 
The acceptable fit range is defined as being compatible with at least all but one of the 95% C.L. regions in fig ||. The 
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FIG. 3: Lepton Et distribution at the best-fit point in the data (points), SM background (dashed histogram) and our best-fit 
point (solid histogram). \/iV uncertainties have been imposed upon the data. 



relevant branching ratios of the smuon are 

BR{i1l Xl^J■) = 0.984, BR{J1l ~^ ud) = 0.015, BR{ilL i^G) = 0.001, (2) 

with a Ufctimc of 1 x 10^^^ sec, whereas for the lightest neutralino we have 

BR{xi G-f) = 0.975, BR{xi ^ Ge-e+) = 0.020, (3) 

with a lifetime of 1 x 10^^^ sec. At such small values of \2ii and m^,, R-parity violating decays of the lightest 
neutralino are negligible. In Table 0, we show the percentage of events making it through each of the cuts. The table 



TABLE L MSSM spectrum used to explain anomalous events at the best fit point within the acceptable fit range for mg, = 10"'^ 
eV, tan/3 = 10 and A211 ~ 0.01. We have displayed the relevant sparticle masses. All other sparticles are heavier than 1900 
GeV, and thus do not interfere with our analysis. 

particle 



best-fit mass 
range 



X? 



IJ.R,eR 



131 GeV 104 GeV 87 GeV 130 GeV 
121-162 GeV 92-141 GeV 87-120 GeV 120-161 GeV 
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TABLE II: Percentage of SUSY events for the best fit point that satisfy cumulative cuts for /i7-^y events at CDF, Run I. 
Events that pass a cut in a given entry also pass those cuts to the left. 



cut 


Eri-f) > 25 GeV 


isolated detected 7 


£t(m) > 25 GeV 


M < 1-0 


detected /i 


^j. > 25 GeV 


percentage 


80.8 


48.6 


35.2 


22.8 


13.3 


10.9 



shows that 11.4% of the smuons produced end up as detected fJ,j^T events in CDF. The corrected cross-section of 
0.091 pb corresponds to 7.8 events additional to the 4.2±0.5 predicted by the SM for 86 pb^^ of luminosity, adequately 
fitting the excess of events quoted by CDF at Run I. 

We now determine the rate of single sneutrino production at Run I. The process is: v —> vxi followed by Xi Gj. 
This would appear to mimic Zj production, where Z vD. To compute the cross-section for this process, we use the 
cuts used by the DO experiment in their 7^7- analysis |l^. With their cuts, we predict a supersymmetric cross-section 
of 0.054 pb for the .^7-7 process at the Run I energy, which corresponds to about 0.7 events for the 14 pb"^ data 
analyzed by the DO experiment. The DO experiment observed 4 events over a SM background of 1.8±0.2 events but 
with a much bigger background coming from cosmic ray sources which is estimated to be 5.8±1.0. As far as we are 
aware, the analysis has not yet been done with the full Run I Tcvatron data but we would expect about 5.4 events 
for a 100 pb~^ data sample. 

We perform the above analyses for Run II (at ^/i = 2 TeV) for the best fit point in order to make predictions for 
observable supersymmetric cross sections: 

aiju^j.) = 0.098 pb, (jh^r) = 0-36 pb, (4) 

which, ought to be observable with good statistics. Since the numbers for the cuts and the efficiencies at Run II 
are not available, we have simply used those that the CDF experiment used in their /i7-^y analysis at Run I. To 
that extent, these numbers are only indicative. For example, with an integrated luminosity of 2 fb^^, these cross- 
sections would correspond to 195 and 720 events, respectively. We predict 0.8 expected selectron pairs at Run I. 
Thus, the discrepancy with respect to the SM ||^ from the observation of an eejj^rp event in the Run I data is 
vastly ameliorated. R-parity conserving production processes such as these will be observable at Run II providing 
more independent checks upon our scenario. One expects an identical number of smuon pairs, leading to a fifi'yj^j' 
final-state. This final state has not yet been observed by CDF, but we note that combining the eejj^rp and fi^jj^rp 
channels, our model still vastly ameliorates the discrepancy with respect to the SM. 



III. CONCLUSIONS 



We have demonstrated that R-parity violating supersymmetry with a light gravitino can explain an anomalously 
high measured cross-section for the fJ,j^T channel. We have provided possible tests for this hypothesis, in the form 
of SUSY cross-sections for the 7-^^ channel, and predictions for the cross-sections of both channels at Run II of the 
Tevatron collider. The 7-^^ channel looks particularly promising because it will allow an independent check of our 
scenario. 

Another interesting question to ask is whether the signal can also be obtained from a specific model of supersym- 
metry breaking consistent with all other data. 
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